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, equation of the time-dependent effective Dirac Hamiltonian is studied in de- 
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, time, a relation is presented between the spectral flow and two characteris- 



tics of the background gauge field. These characteristics are the well-known 



"winding factor," which is defined to be the change of the Chern-Simons num- 



ber of the associated vacuum sector of the background gauge field, and a new 
"twist factor," which can be obtained from the zero-eigenvalue equation of 



H ' the effective Dirac Hamiltonian but is entirely determined by the background 

gauge field. For a particular class of (weakly dissipative) Liischer-Schechter 
gauge field solutions, the level crossings are calculated directly and nontrivial 
contributions to the spectral flow from both the winding factor and the twist 
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to electroweak baryon number violation in the early universe. 
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I. INTRODUCTION 



An important consequence of the triangle anomaly |I|J^ is the violation of fermion number 
conservation in the electroweak Standard Model if the background gauge field has nonva- 
nishing topological charge Q. This connection between the triangle anomaly and fermion 
number violation in the electroweak Standard Model was first pointed out in Ref. 0. For 
the background gauge field, the calculation of Ref. used Euclidean instanton solutions 
in order to calculate the tunneling amplitude between topologically different vacua. A 
special feature of the Euclidean (imaginary-time) approach is that gauge field configurations 
with finite action fall into homotopy classes labeled by an integer Q. This integer Q then 
gives the number of fermions produced. 

For real-time fermion- number- violating processes sphaleron-like gauge field config- 
urations are believed to play crucial role p. These gauge fields "interpolate" between 
topologically different vacua and have sufficiently high energy to overcome the energy bar- 
rier. But, in contrast with the Euclidean approach, the anomalous fermion production in 
Minkowski spacetime is not, in general, given by the topological charge Q of the classical 
gauge field background. The reason is that Q may be a noninteger or even, for the case of 
Yang-Mills-Higgs theory, not well-defined; cf. Refs. It is not clear which quantity, in 

general, determines the anomalous fermion production for real-time processes. 

For pure Yang-Mills theory in Minkowski spacetime, the authors of Refs. P,p!0| have 
argued that the number of produced fermions is given by the change of winding number of the 
associated vacua of the initial and final gauge field configurations. The "associated vacua" 
of a given classical gauge field background represent the particular vacuum configurations 
that the background field would approach as t — *• — oo or t — >■ -|-cxd. As is clear from the 
context, the authors of Refs. P,p!0| considered dissipative background fields in order to be 
able to quantize the fermion fields at t = ±cxd. 

The question, now, is what happens to fermion number violation if a powerful energy 
source creates a nontrivial gauge field background over an extended spacetime region (for 
example, in a high-energy collision experiment or in the early universe). In this case, we 
cannot readily associate the initial high-energy state with a particular vacuum configuration 
and the field-theoretic approach used in Refs. breaks down. For nontrivial classical 

bosonic background fields, it is, moreover, not known how to construct a fermion number 
operator in terms of the quantized fermionic fields . 

Still, fermion number violation can be directly observed from the level crossing of the 
energy eigenvalues of the time-dependent Dirac Hamiltonian, see Ref. and references 
therein. The overall effect of level crossing can be characterized by the "spectral flow," 
defined to be the number of eigenvalues of the Dirac Hamiltonian that cross zero from below 
minus the number of eigenvalues that cross zero from above, for a given time interval and 
direction of time. 

In this paper, we study the zero-eigenvalue equation of the effective Dirac Hamiltonian 
for spherically symmetric chiral fermion fields and classical SU(2) Yang-Mills gauge field 
backgrounds. A relation is found between the spectral flow and certain features of the 
spherically symmetric background gauge field. These features are the well-known topological 
"winding factor" and a new type of "twist factor," both of which will be defined later. 
The spherically symmetric subspace of (3+l)-dimensional chiral SU{2) Yang-Mills theory 
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is equivalent to an (l+l)-dimensional U{1) gauge theory coupled to a Higgs-like complex 
scalar field ]ll]| and several two- component Dirac fields. This drastic simplification allows 



us to examine the problem using analytical methods. 

The (3-|-l)-dimensional spherically symmetric SU{2) gauge field backgrounds considered 
in this paper are, in general, nondissipative, which means that the energy density does not 
approach zero uniformly as t ^ ±00. For these nondissipative gauge field backgrounds, 
the contribution of the twist factor to the spectral flow is manifest. The general result for 
the spectral flow also applies to dissipative spherically symmetric gauge field backgrounds, 
which will be classified later as "weakly dissipative" and "strongly dissipative." It will be 
shown that the nonvanishing effect of the twist factor for the spectral flow can already 
appear for weakly dissipative spherically symmetric gauge field backgrounds, such as certain 
Liischer-Schechter gauge field solutions [0,0. Our paper may, therefore, be viewed as a 
continuation of the work of Refs. ||5|,p|,p!0| 



The outline of this paper is as follows. In Section II, we present the model and the basic 
formalism. After giving the chiral SU{2) Yang-Mills theory in the spherically symmetric 
Ansatz, we briefly review the topological properties of the gauge field background. The 
gauge field topology is, in the first place, characterized by the winding factor, defined to be 
the change of the Chern-Simons number of the associated vacuum sector of the gauge field 
configuration. 

In Section III, we consider the zero-eigenvalue equation of the time-dependent effective 
Dirac Hamiltonian. By investigating the zero-eigenvalue equation directly, we are able to 



identify a family of Riccati equations ||1J,|T5[, from which the twist factor of the spherically 
symmetric SU (2) gauge field configuration can be obtained. 

In Section IV, we present a result for generic spherically symmetric SU{2) gauge field 
backgrounds, which relates the spectral flow to both the winding factor and the twist factor. 

In Section V, we investigate the level crossing phenomenon for the particular spheri- 
cally symmetric SU{2) gauge field backgrounds given by certain Liischer-Schechter solutions 
||l^ , p!3[| and verify our relation for the spectral flow. Specifically, we demonstrate the signifi- 
cant effect of the twist factor for a class of Liischer-Schechter background gauge fields with 
energies far above a sphaleron-like barrier. 

In Section VI, finally, we summarize our results and briefly discuss the role of dissipation. 
There is also an appendix, which provides the proof of a result needed in this section. 

For the benefit of the reader, we remark that Sections II B, III C, and IV B form the 
core of the paper. 



II. CHIRAL SU{2) YANG MILLS THEORY 

In this section, we review the spherically symmetric Ansatz for massless chiral fermions 
coupled to classical SU{2) Yang-Mills gauge fields and establish our notation. Furthermore, 
we recall the definition of the topological winding factor. 



A. Spherically symmetric Ansatz 
The SU (2) Yang-Mills theory with massless chiral fermions is described by the action. 
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S = Sg + S, 



^G = -lhl d^xTliF^'^Frnn). Sf= [ d^X ^^V]/^, (2.1) 



where Sq represents the gauge field action and Sp the fermionic action. Latin indices m, 
n, etc. run over the coordinate labels 0, 1, 2, 3, and the metric tensor for fiat Minkowski 
spacetime is rjmn = diag(— 1, 1, 1, 1). Repeated indices are summed over. The fiavor index 
/, in particular, is summed over 1, . . . , Np. Also, natural units are used for which c = h = 
1. 

The SU{2) field strength tensor Fmn and the covariant derivative Dm for the fermionic 
fields are defined as follows: 

Fmn = dmA-n — OnAm + [Am, An], Am = A,^ T / (2i) , 

Dm = dm + AmPL, Pl = (1 " Ts) /2, = (1 + Tg) /2 . (2.2) 

The Dirac matrices are taken in the chiral (Weyl) representation, 

-(J?)' ^'-(li)' ^^-(°oO' "'-{l-.)- 

The conjugate spinor is given hy = \l/j(— iF^). Here and in the following, r" and cr" 
are Pauli matrices carrying isospin and spin indices, respectively. The action ( p.l|) thus 
corresponds to a chiral SU{2) gauge theory, with interacting left-handed fermions {"^if = 
Pl^ f) and noninteracting right-handed fermions rj = Pr^ f). 

The total number Np of fiavors in the fermionic action ( |2.1| ) must be even, in order to 
cancel the nonperturbative SU{2) anomaly [Q. Henceforth, we focus on a single fiavor and 
drop the index /. Since there is no natural mass scale for the classical SU{2) Yang-Mills 
theory, we also take an arbitrary mass scale to work with. (In the full theory, quantum 
effects may, of course, fix the scale.) 

In this paper, we concentrate on the spherically symmetric subspace of the (3+1)- 
dimensional theory. We use the following Ansatz for the gauge fields: 

Ao(x) = ^ao(t,r)r-x, (2.4a) 



A.(x) = 1 



ait, r) — 1 „ p(t,r) „ „ ^ , . ,^ 
eabcXcU H [dab - XaXb)n + ai{t,r) (r ■ X)Xa 



(2.4b) 



where a, (3, oq, and ai are real functions of t and r; cf. Refs. [pT] , p^ . These gauge fields are 
invariant under spatial SO (3) rotations, up to a gauge transformation 

Am^AZ = n{dm + Am)n\ (2.5) 

with n{x) G SU{2). 

The spherically symmetric Ansatz for the fermionic fields is given by (see Ref. [l^ and 
references therein) 
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r2 ^^(x) 



iGL{t,r)T- x] 



-1 

Wo' 



V2 



■ X 



v 






+1 



isospin 



isospm / spin 



isospm 



(2.6a) 



(2.6b) 



:2.6c) 



isospm / spin 



where H^, Hji, Gl, and G/j are complex functions of t and r. In components (a for isospin 
and a for spin), the two constant spinors of Eqs. (^]6|b,c) can be written as eaa and 6aa, 
respectively, where e and S are the Levi-Civita and Kronecker symbols. 

Furthermore, we assume that all physical (3+l)-dimensional field configurations are de- 
scribed by infinitely differentiable functions (this assumption can be relaxed). In order to 
have regular behavior at the spatial origin for the (3+l)-dimensional field configurations and 
their derivatives, the Ansatz functions should satisfy the following r-parity expansions near 
r = 0: 



ao(t,r) 
a(t, r) 
HL,R{t^r) 



2k+l 



1 + E 



k=i 



a 



(2k) 



{t)r 



2k 



TT(2fc) 



{t)r 



2k 



ai(t,r) 
/3(t,r) 
G'L,fi(t,r) 



{2k+l)i^-^ _^2k+l 



k=0 '-^L,R 



(2.7) 



with the expansion coefficients depending on time only. 

If we substitute the Ansdtze (p.4| ) and (|2.6|) into the action ( p.l|) , the following reduced 
actions are obtained (for a single fermion flavor): 



G 



Sf 



An 
An 



dt 



dt 



dr 



r'Uur + \D,x\' + j^,{\x\'-l) 



(2.8a) 



dr <^ YD^ + - (Rex + ^75 Imx) + I'^d^ + 



(2.8b) 



Greek indices /i, u, etc. run over the coordinate labels 0, 1, and are lowered with the metric 
rj^u = diag(— 1,1). The coordinates (x°,x^) correspond to (t, r). The theory ( |2.8D can be 
interpreted as an (l+l)-dimensional U{1) gauge field theory with a Higgs-like complex scalar 
field xitjf) and two-component Dirac spinors \l/i(t,r) and \l/r(t, r). In terms of the Ansatz 
functions, the U{1) field strength f^^^, the complex scalar and Dirac fields, and the covariant 
derivatives are given by: 
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^,,,^v[/t^(-^7°), D,^i^{d, + z{a,/2)j,)^i, (2.9) 

with 

70 = za\ 7' = -^^', 75 = -7V = cr^ (2.10) 



The spherically symmetric Ansatz ( |2.4| ), ( |2.6D preserves a f/(l) subgroup of the SU{2) 
gauge group, with transformation parameters i7(x, t) = exp[iuj{t,r) r ■ x/2] in Eq. ( p.5| ) 
above. Under these particular SU{2) gauge transformations, we have for the (1+1)- 
dimensional fields the following U{1) gauge transformations: 



a 



a^ + d^u, x^e^'x, ^ e-^("/2)^^ ^z, ^ (2.11) 



In order to maintain the regularity of the (3+l)-dimensional field configurations, uj{t,r) 
should have an odd r-parity expansion near r = 0, 



oo 

u;{t,r) = Y,^^"'^'\t)r''+\ (2.12) 

A;=0 



where the expansion coefficients are dependent only on time. 

For later reference, the (l+l)-dimensional fields with finite energy approach a vacuum 
configuration at infinity, provided 

X^e^", a^^d^u, ^ 0, ^z,. ^ 0, for r ->cx) , (2.13) 

and are regular at the spatial origin, provided 

IxHl, D^x^O, ^i,.^0, for r^O. (2.14) 



See Ref. |jT^ for further details. Throughout this paper, we consider regular spherically 
symmetric SU (2) gauge fields with finite energy. 



B. Gauge field winding factor 

For the description of the topology of spherically symmetric SU{2) gauge field back- 
grounds, it is convenient to express the (l+l)-dimensional complex field x(^)^) iii polar 
form: 

X{t,r)=p{t,r)exp[iy^{t,r)], p(t, r) > 0. (2.15) 

The "associated vacuum sector" of the background gauge field at a fixed time t is obtained 
from the configuration with p{t,r) replaced by 1, but with (p{t,r) and a^{t,r) unchanged. 
[Note that the resulting configuration with p(t, r) = 1 may still have nonzero energy density 
( p.8a| ).] For the gauge choice x(^) 0) = oo) = 1, the integer winding number is then 
defined as 
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N^{t) = [v{t, oo) - v{t, 0) ]/(27r). (2.16) 

This winding number N^{t) is, in fact, equal to the Chern-Simons number of the p = 1 
gauge field at time t; see Eq. (2.18) of Ref. |T^. For a particular x(^;''^) configuration, the 
winding number N^it) is unambiguous, provided |x(t,r)| > 0. See Sec. IV B for further 
discussion. 

For a time interval [ti,tf] with ti < tf, generic spherically symmetric SU{2) gauge field 
backgrounds are characterized by the change of winding number A^^ between the initial and 
final configurations, 

AN^[tf,U]^N^{tf)-N^{U). (2.17) 



Henceforth, we call AN^[tf,ti] as defined by Eq. ( |2.17| ) the "winding factor" of the spheri 



cally symmetric SU{2) gauge field. Our definition of the "winding factor" is directly inspired 
by the results of Ref. [^], obtained for a particular class of background fields that will be 
discussed further in Sec. V. 



III. TIME-DEPENDENT DIRAC HAMILTONIAN AND TWIST FACTOR 

In this section, we consider the zero-eigenvalue equation of the time-dependent effective 
Dirac Hamiltonian for a given spherically symmetric SU{2) gauge field configuration at one 
particular time. The existence of fermion zero modes is discussed and a necessary condition 
derived. In addition, the so-called twist factor is introduced, which will play an important 
role in Sec. IV. 



A. Fermion zero modes and level crossings 



The general solution \l/;(t,r) of the (l-l-l)-dimensional Dirac equation from the action 
( p.8b|) can be expressed as a linear combination of the eigenfunctions of the corresponding 
time-dependent Dirac Hamiltonian. The eigenvalue equation of this Hamilton operator is 



H{t,r) ^(t,r) = E{t) ^(t,r), 

H{t, r) = 75 ao/2 - ij^ Di 

+ ^7° (Rex + i Imx 75) / t. 



(3.1a) 
(3.1b) 

where the covariant derivative Di has been defined in Eq. (|2.9| ) and \I' now stands for 
the two-component Dirac spinor of that same equation. [The other Dirac field \l/r(t,r) 
of the action ( |2.8b|) has no interactions and will not be considered in the following.] The 
Hamiltonian ( |3.1bD depends on t and r through the background fields xit^''^) '^fii^i'^)) 
together with an explicit dependence on r in the Z7° term. 

It is known that the zero-crossing of an energy eigenvalue of the Dirac Hamiltonian is 
one of the crucial ingredients of fermion number violation; cf. Refs. In our case, the 

zero-eigenvalue equation (|3.1a]) at fixed time t can be written as 
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dr'^=A'^, A = Ah + Aa, (3.2a) 
Ah = - (7^ Rex + i7° Imx) /r, (3.2b) 
^^ = -200/2-27501/2, (3.2c) 

where 9^ stands for the partial derivative with respect to r. For later convenience, we have 
decomposed the matrix A of Eq. (|3.2aD into a Hermitian part Ah and an anti-Hermitian 



part Aa- Recall that we use two-dimensional Dirac matrices 7° = ia^, 7^ = —cr^, and 
75 = 0"^, with 0"" the standard 2x2 Pauli matrices. 

In order to have a regular (3+l)-dimensional fermionic field at r = 0, the (1+1)- 
dimensional fermionic field \E'(t, r) must satisfy the boundary condition ^^(^,0) = 0, which 
is already implemented by the Ansatz ( |2.9| ). A fermion zero mode is then defined to be a 
normalizable solution of Eq. (|3.2a|) with boundary condition \E'(t,0) = 0. Specifically, the 
normalization condition is given by 

dr |^(t,r)|2 = 1. (3.3) 

The existence of a fermion zero mode at a particular time t = t* does not necessarily 
imply level crossing of the eigenvalue of the Dirac Hamilton operator. In fact, the energy 
level E{t) could just "touch" the E = Q value instead of "crossing" it. Therefore, it is 
necessary to check that level crossing really occurs. This can be done by calculating the 
time-gradient of the energy eigenvalue E{t) at t = t* . If (\E/(lt\t=t* 7^ 0, then there is level 
crossing ai t = t* . 

The overall effect of level crossings can be characterized by the "spectral flow" 
defined to be the number of eigenvalues of the Dirac Hamiltonian H{t,r) that cross zero 
from below minus the number of eigenvalues that cross zero from above, for the time interval 
[ti,tf] considered. The spectral flow will be discussed further in Sec. IV. Here, we continue 
the investigation of the zero-eigenvalue equation per se. 



B. Gauge-invariant zero-eigenvalue equation 



We first express the zero-eigenvalue equation ( p.2aD at a fixed time t in terms of a set of 
bosonic background fields that are invariant under the U{1) gauge transformations ( p. 111) . 
As can be seen from Eq. (|3.2|) , a nonvanishing gauge field ao contributes only a complex 
phase factor to the solution \I'. Up to an overall phase factor, the solution is then 



\E'(t, r) = exp 



i / dr' ao{t,r')/2 



ro 



^M(t,r) 



provided "^uityr) satisfies the nontrivial matrix equation 

dr^M =[-il5 ai/2 - (7^ Rex + «7° Imx) /r ] ^ 
with boundary condition 

^M(t,0) = 0. 



M, 



(3.4) 



(3.5) 



(3.6) 



8 



The existence of a fermion zero mode for the hnear differential equation ( p.2|) is thus equiv- 
alent to having a normalizable solution of Eq. (|3.5| ) with boundary condition (|3.6|) . 
Next, we apply a unitary transformation to Eq. ( p.5| ), 

^M^^A^A^^M, (3.7) 

with the transformation matrix 

A = -2 7°exp[iy.75/2], (3.8) 

that diagonalizes the Hermitian matrix Ah via Ah — * A^Ah-^- We obtain the following 
zero-eigenvalue equation for \1/a: 

9,*A = (A + A)^A, (3.9a) 

Ao = A^AhA = \i^ = -\ a^ (3.9b) 

Ai = A^ {-dr - A = m-if, =7^^c^^ (3.9c) 

with the further definitions 

A = p/r > 0, 7^ = (ai - dr^)/2. (3.10) 



It follows immediately from the definition ( p.l5| ) that the matrices and Ai are invariant 



under the U{1) gauge transformations ( p.ll| ). Moreover, and Ai are real matrices and 



the solution \1/a can be taken real, up to an overall complex phase factor. In the following, 
we take \I/a to be strictly real and drop the subscript A. 

For finite-energy background gauge fields with r-parity expansions as given by Eq. ( p.TD , 
one can show that the following limits hold: 

lim7^/A = lim 7^/A = 0. (3.11) 

r— >0 r— ►oo 

This demonstrates that the diagonal matrix determines the local structure of the solution 
of the differential equation (|3.9| ) in the regions of small and large r (see also Ref. [1T9| ). 



C. Spinor twist number and twist factor 

Since the solution of the transformed zero-eigenvalue equation ( p. 91 ) is taken to be real, 
one can write ^ in polar notation, 

m{t,r) = \m{t,r)\ exp[z75e(t,r)] (^J^ , (3.12) 

where G M measures the relative rotation of the spinor away from the \I^2-axis in the 
configuration space of Recall that 75 = cr^, so that the exponential factor in Eq. ( |3.12| ) 
reads U cos + ia"^ sin 0. 

From Eqs. ( |3.9| ) and ( |3.12| ), one finds that and at fixed time t satisfy the following 
coupled differential equations: 
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with the definitions 



drO = P[e] +7^, 
ai^l = A 1^1 COS 29, 



V[e] = -X sin 29, 7^ = (ai - 9^<^)/2, 
A = p/r > . 



(3.13a) 
(3.13b) 



(3.14) 



In order to obtain regular behavior at r = 0, the solutions of the differential equations (|3.13|) 
must satisfy the following boundary conditions: 



9(t,0) = mod TT, 
1^(^,0)1 = 0. 



(3.15a) 
(3.15b) 



More specifically, these boundary conditions are needed because A is singular at the spatial 
origin r = 0; see Eq. ( ^.141) . 

The nonlinearity of the differential equation ( ^.13a| ) originates from the fact that the 
linear differential equation (|3.9|) mixes the components of the spinor Furthermore, the 
differential equation ( p.l3a| ) involves only 9, whereas Eq. ( ^.13b| ) contains both 9 and |\&|. 
These two properties of Eq. ( p.l3a| ) will turn out to be crucial for the results of the present 
paper. 

Remarkably, the nonlinear differential equation ( |3.13a| ) for a given time slice t can be 
transformed into a generalized Riccati equation |1T^JT5[| by setting Y{t,r) = tan9(t,r), 

drY -n{l + Y^) +2XY = . (3.16) 

The analysis is, however, best carried out with the nonlinear differential equation in the 
form as given by Eq. ( |3.13aD , where the term V is called the "deviator" and the term 
TZ the "rotator," for reasons that will become clear shortly. Henceforth, we refer to the 
single differential equation ( p. 13a ), with the implicit boundary condition ( p. 15a ), as the 
"transformed Riccati equation." 

Let us consider the asymptotic behavior of the solution Q{t, r) of the transformed Riccati 
equation ( p.l3a| ) at a fixed time t. The deviator T> dominates, in general, the right-hand side 
of Eq. (|3.13a|) for large r, according to Eq. (|3.11|) . For large r, Eq. (|3.13a|) can therefore be 
approximated by 



dr& = -A sin 29. 

The differential equation (|3.17|) has three types of solutions at a fixed time slice t, 

e{t,r) = Ntt, 
9(t,r) = (iV'+ l/2)7r, 
tan[9(t,r)] = tan[9(t,ro)] 



X exp 



-2 



dr'A(t,r') 



ro 



(3.17) 

(3.18a) 
(3.18b) 

(3.18c) 
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for arbitrary integers and A^'. 

The nontrivial solution B(t, r) given by Eq. ( 3.18c| ) is attracted toward the value A^vr as 
r oo, since A(t, r) is non-negative and has a divergent integral toward infinity. This shows 
that the "point" Q{t,r) = Nir, with G Z, is asymptotically stable in the solution space of 
the differential equation ( pTzD ; cf. Ref. 191. For the trivial solution Q{t,r) = {N' + 1/2) tt, 
an arbitrarily small deviation will lead to a nontrivial solution given by Eq. ( p.l8c|) , which 
asymptotically approaches the value N" tt, with A^" G Z. The "point" G(t,r) = (A^' + 
1/2) TT, with A^' G Z, is thus asymptotically unstable in the solution space of the differential 
equation (|3.17|). 



The solutions of the complete differential equation (|3.13a| ) with boundary condition 
( p.l5a| ) can thus be classified according to their asymptotic behavior. At a fixed time slice 
t, there are two classes: 



SN{t) = {Q{t,r)\ lim e(t,r) = A^vr}, 

r— ►oo 

UN'it) = {e(t,r)| lim e(t,r) = (A^' + 1/2) tt}, 



(3.19a) 
(3.19b) 



with N,N' G Z. If 6 G Sn (© G Un), then 9 is asymptotically stable (unstable) in the 
solution space of the differential equation ( p. 13a] ). 

For any solution B(t,r) G UN{t) with arbitrary integer A^, there necessarily exists a 
fermion zero mode, as follows from Eqs. ( 3.13b| ) and ( |3.15b|) . [It is clear that the normal- 
izability condition (|3.3|) of the fermion zero mode requires the asymptotics of Eq. (|3.19b| ).] 
Therefore, it suffices to study the transformed Riccati equation ( p.l3a| ) in order to determine 
the existence of a fermion zero mode at a particular time t. 

At this moment, we can explain the use of the terms "deviator" and "rotator" in the 
transformed Riccati equation ( 3.13a|) . The observation from Eqs. ( 3.17]) and ( p.l8c ) is that 
V pulls B(t,r) toward the value A^vr as r — oo. In other words, it leads to a deviation 
of G(t,r) from the special path approaching the value (A^' + 1/2) tt as r — > oo, for which 
a fermion zero mode exists. This is the reason for calling the term T> in Eq. (|3.13a|) the 
"deviator." In the absence of the deviator V over the interval [vq, ri], say, one observes from 
Eq. ( |3.13a ) that TZ generates a simple rotation of the spinor by the angle AG = /J'' dr'7^(r'). 
This is then the reason for calling the term TZ in Eq. p.l3a| ) the "rotator." The fermion 
zero mode solutions will be discussed further in the next subsection. Here, we continue the 
discussion of the transformed Riccati equation from a more general viewpoint. 

Using analyticity and the Cauchy-Lipschitz existence and uniqueness theorem for or- 
dinary differential equations [14,|T5il9|, it can be shown that the solution B(t,r) of the 
transformed Riccati equation ( p. 13a] ) with boundary condition O(t,0) = is unique. For 
the regular finite-energy gauge fields considered, it can also be shown that the solution 
Q{t, r) is bounded. 

The uniqueness of the solution Q{t,r) and its asymptotic behavior allow us to classify 
the gauge field background at one particular time t by the quantity 



Ne{t) = [Q{t,oo)-Q{t,0)]/7T, 



(3.20) 



which can take integer or half-odd- integer values. According to the definition ( ]3.12] ), the 
mapping G(t,r) = exp[z 75 9(t, r)] G 5*0(2) gives the twisting of the spinor in the configu- 
ration space of \I' for fixed time t. We therefore call NQ{t) the "spinor twist number." 
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It is now convenient to characterize a time- dependent spherically symmetric SU{2) gauge 
field background by the change of spinor twist number between initial and final configura- 
tions. Henceforth, we call 



ANe[tf,U] = Ne{tf)-Ne{U 



(3.21) 



the "twist factor" of the spherically symmetric SU{2) gauge field over the time interval 
[ti,tf], with ti < tf. 

It is important to realize that the twist factor AA'^e measures an intrinsic property of 
the SU{2) gauge field configuration. Formally, Eqs. (|3.9| a), (|3.12| ), ( p.20| ), and (|3.21| ) give 



ANe[tf,ti 



TT 



dr 



dt 



d ( d 



dt V dr 



-Q{t,r] 



Q(t,r) = — Tr I ia"^ In < lim ePexp 



dr' {A(t,r')+A(t,r')} 



(3.22a) 



(3.22b) 



where V represents path ordering. Here, ^0(^5''^) and Ai{t,r) are defined by Eqs. (p.9|b,c), 
in terms of the (l+l)-dimensional gauge field functions p{t, r), (p{t, r), and ai(t, r). Whether 
or not there exists a more direct way to obtain AA'^e remains an open question. 



D. Necessary condition for fermion zero modes 

With the results of the previous subsection, it is possible to find a necessary condition for 
the existence of fermion zero modes at a particular time t. We first introduce the following 
diagnostic: 

/C±(t) = [ dr7^(t,r) 

POO 

= dr^[±7^(^,r)] 7^(^,r) , (3.23) 
Jo 

with the domains of positive or negative values of TZ{t, r) defined by 

D±(t) = {r| sgn[7^(t, r)] = ±1} C [0, 00) (3.24) 

and 9 the usual step function, 9[x] = for x < and 9[x] = 1 for x > 0. Note that, by 
definition, /C+ > and /C_ < 0. Note also that the rotator 7l(t, r) from Eq. (|3.14| ) is entirely 
defined in terms of the background fields ai(t,r) and (p(t,r) = arg x(t,'>^)- 

Consider the transformed Riccati equation (p. 13a]) with boundary condition Q{t, 0) = 0. 
The integration of TZ{t,r) over the domain D^{D^) then accounts for the rotation of the 
spinor in the "+"("—") direction. But the deviator V, for values B G (— 7r/2, 7r/2), brakes 
the rotation forced by the rotator TZ. The crucial point, now, is that in order to have a 
fermion zero mode at time t the total action of the rotator should overcome the resistance 
from the deviator in the region —tt/2 < O < tt/2, so that the solution G(t,r) ends up with 
|0|>7r/2atr = 00. A necessary condition for the existence of a fermion zero mode at one 
particular time t is, therefore. 
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/C„,ax(t) = max [/C+(t) , I /C„(t) I ] > tt/2. (3.25) 



Having established the necessary condition ( p. 25 ), it would certainly be interesting to 



obtain also a necessary and sufficient condition for the existence of a fermion zero mode in a 
given static gauge field background. But, without further input, it appears difficult to find 
such a condition. For this reason, we turn in the next section to the role of time-dependent, 
continuous gauge field backgrounds. 



IV. SPECTRAL FLOW 



In this section, we consider the spectrum of the effective Dirac Hamiltonian ( |3.1b| ) for 



time-dependent spherically symmetric SU{2) gauge fields. In Sec. IV A, we derive a relation, 
Eq. ( [4. 321 ), between level crossing and the change of winding number or spinor twist number 
over an infinitesimal time interval. From this result, we obtain in Sec. IV B the appropriate 
relation, Eq. ( [4.39|) , for the spectral flow over a finite time interval. Section IV A is rather 



technical and may be skipped on a ffist reading. 



A. Level crossing from changes in winding and twist numbers 



1. Perturbative expansion 



We start from the transformed Riccati equation ( p.l3a| ), with boundary condition ( |3.15aD 



at a particular time t = t* and study the change of the solution 0(t, r) in the neighborhood 
of t = t*. For finite r and t = t* ± e with e an arbitrarily small positive constant, one can 
expand the background fields as follows 

X{t* ± e, r) = A(r, r) ± e 9tA|,^,^ + 0{e^), (4.1a) 
7^(r ± e, r) = 7^(r , r) ± e 9^7^|,^,^ + 0(6^), (4.1b) 

where the upper and lower time derivatives of the background fields are defined by 

9tA(t,r)|^^^^ = Jim 9tA(t,r), 

dtKt.r)\,^t- =limatA(t,r), (4.2) 

and similarly for dtTZ. The solution 0(t, r) for t = t* ± e can be written as 

e(r ± e, r) = e±(r , r) ± e h{t\r) + 0{e^). (4.3) 

The function /i(t, r) is continuous at t = t*, but the ffist term on the right-hand side of Eq. 
( [4. 3D allows for a discontinuity. For the moment, we consider the functions 0±(t*,r) in Eq. 
( [4.3| ) to be equal. 

By substituting Eqs. ( [4.1D and (|4^ ) into the transformed Riccati equation (|3.13a|) , one 
obtains to ffist order in e the following linear differential equation for fi\ 
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[dr + 2 A(f , r) cos 2e(r , r) ] /i(t* , r) = r), 

with the definition 



(4.4) 



ji(r,r) = 9^7^|,^,, - 9tA|,^,, sin 29, 

and boundary condition 

/i(r,o) = o. 

The solution of the differential equation (|4.4| ) is found to be 

h{t\T)=J,{t\T)/\^{t%T)\\ 

with the definition 

Ji{t\r)= r dr'\^it*,r')\^nit*,r'). 
Jo 

Here, we have used the solution of Eq. ( p.l3b|) . 



l^^it, r) \ oc exp 



dr A(t,r ) cos20(t,r') 



L ^ ro 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



to obtain Eq. ( [4.7|) in the form shown. 

The function ji{t,r) is continuous at t = t* for smooth background fields X{t,r) and 
TZ{t,r), which implies that the solution /i(t, r) of Eq. ([4.4|) is also continuous at t = t*. If, 
on the other hand, the time slice t = t* corresponds to a local change of the gauge field 
winding number A^^, the partial derivatives of X(t,r) and 7l(t,r) are not well-defined at 
t = t*, as will be shown later. In this case, the function ji{t,r) is not well-defined either, 
which affects the continuity of the solution fi{t,r) of Eq. ( |4.4| ). In Sec. IV A 3, we will 
show that the function f\{t,r) can be taken to be continuous at t = t*, provided possible 
discontinuities are accounted for by the leading terms G_(t*,r) and 0+(t*,r) in Eq. ( |4.3|) . 



2. Time-differentiable A and TZ 

Consider a particular time slice t = t*, for which the A(t,r) and TZ{t,r) fields are differ- 
entiable with respect to time, 

dt\it,r)l^,,^ =dt\it,r)\,^,,_, (4.10a) 

a,7^(t,r)|,_ =a,7^(^,r)|,^,. , (4.iob) 



with upper and lower time derivatives as defined in Eq. 

First, suppose that there is no fermion zero mode at t = t*, so that <d(t*,r) ^ N tt as 
r —>■ oo. For large r and using Eq. ( p. 11 ), one then obtains from Eq. ( |4.4| ) the results 
lim^^oo ji/A = and /i oc r~^, which imply 

lim /i(r,r) = 0. (4.11) 
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This indicates that there is no change of the asymptotic behavior of Q{t,r) in the neigh- 
borhood of t = t*, which corresponds to having a constant spinor twist number ( p. 20 ) at 
t = t\ 



e\t=t' 



= Ne{t* + e) - Ne{t* - e) 



0. 



(4.12) 



[We reserve the notation AA^^e for the global change of spinor twist number; see Eq. ( ^■21[ ). 
Of course, ^A^^e is in no way "infinitesimal," see Eqs. ( |4.14|) below.] 

Next, consider the case of having a normalized fermion zero mode at t = t* with 
Q{t*,r) E Un, that is, belonging to the "unstable" class of solutions ( p.l9b| ). Since the so- 
lution Q[t*,r) e Un is asymptotically unstable, one observes from Eq. ( |3.18c|) that a small 
positive [negative] perturbation of 6 at large r leads to Q(t*,r) G Sn+i [Q{t*,r) G Sn]- 
From Eqs. (|4.3| ), (|4.7| ), and (^4.8|) one deduces that the fermion zero mode at t = t* sits at a 
bifurcation point for different Nq^s. In fact, the local change of spinor twist number 



SNe\t=t'=Ne{t* + e)-Ne{t*-e) 



is given by 



1 for Ji(r,oo) > 0, 
-1 for Jilt*,oo) < 0, 



(4.13) 



(4.14) 



as long as J7i(t*, oo) ^ 0. [The special case of J7i(t*, oo) = will be discussed in Sec. IV A 
4.] For an elementary discussion of bifurcation theory; see Ref. [ITU . 

Having a fermion zero mode at t = t*, we are especially interested in the time-gradient 
of the fermion energy eigenvalue at t = t*, in order to check for level crossing. The time- 
gradient of the energy eigenvalue of the Dirac Hamiltonian at t = t* is calculated up to the 
first order in e: 



dE 
'dt 



t=t* 



^/(r,r) 
Ji{t*,oo), 



OH 
'dt 



(t* 



(4.15) 



where J7i is defined by Eq. ( [1.8D and '^i{t*,r) represents the (nondegenerate) normalized 
fermion zero mode at t = t* in the two-component spinor notation of Eq. (|2.9| ). The 
expectation value used in Eq. ( [4.15|) is defined by 



{^i{t*,r)\0{t*)\^i{t*,r)) 

POO 

= / dr^iit*,ryOit*)^iit\r) 
Jo 



(4.16) 



for an arbitrary time-dependent Hermitian operator 0{t). 
From Eqs. (I4.14|) and ( [4.15|) , we obtain 



5iVeU. G{-1,+1}, 





'dE 






sgn 


'dt 






with the implicit limit e 


on 


the 


rig 


level crossing and the chan 


ge of 


spinor 



(4.17) 



are time-differentiable at t = t* and generic for times close to it [so that J7i(t*, oo) ^ 0]. 
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3. Time-nondifferentiable A and TZ 



Now, consider a gauge field background for wliicli r) vanislies at tlie spacetime point 
(t*,r*) and tlie winding number A^^ as defined in Eq. ( p.l6| ) changes from a value to 
+ 5N■y^^t=t*■ Generally, the fields A(t,r) and lZ{t,r) are not differentiable with respect to 



time: 



'9t7^(^,r)|,^,. 



(4.18a) 
(4.18b) 



with upper and lower time derivatives as defined in Eq. ( [4.2| ). 

Let us have a closer look at the discontinuities of the time derivatives of A(t,r) and 
TZ{t,r) at t = t*. First, one observes from the definition ( |3.14| ) of the X{t,r) field that the 
time derivative of X(t,r) is not well-defined at the spacetime point (t*,r*): 



= \dtxit*,r*)\/r*. 
Second, introduce the gauge-invariant function R{t*,r) defined by 
R(t*, r) = lim [n{t* + e, r) - 7^(t* - e, r)] 



(4.19) 



(4.20) 



Taylor expanding d^'^it^ r) with respect to both t and r in the vicinity of the spacetime point 
(t*, r*) and using the fact that a{t^ r) — > and /3(t, r) — as (t, r) (t*, r*), one finds that 
R(t*, r) shoots up to infinity at r = r*, whereas it drops to zero for r 7^ r*. Taking the gauge 
condition x(t, 0) = x(^5 00) = 1, one readily proves that 



/ dr R(r,r) = -TT lim^A/"^ 
Jo 

with the local change of winding number defined as 



t=t* 



6N, 



x\t=t* 



N^(t* + e)-N^it*-e). 



(4.21) 



(4.22) 



This shows that the function R(t*,r) is proportional to a Dirac delta-function centered at 
r = r*, 



R(r,r) = -71 lira SN^\t=t* 5{r-r*). 



(4.23) 



The nonvanishing right-hand side of Eq. ( [4.23| ) for r = r* implies that the time derivative 
of the rotator TZ{t,r) is not well-defined at the spacetime point {t*,r*). 

The delta-function-like behavior (|4.23| ) of R can be used to derive the effect of the change 
of the gauge field winding number on the change of the spinor twist number at t = t*. Start 
by defining 



e_(r,r) = lim e(t,r). 

tit* 



(4.24) 
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Then, if 0-(t*, r) belongs to the class Sn, for some integer N, one deduces from Eqs. ( p.l3aD 
O), and (lO^ ) that 



0+(t*,r) = hni e(t,r) 



e_(r,r) -TT lim5iVx|t=t. e(r-r*), 



(4.25) 



where 6 is the usual step function. This implies that the solution changes from one class to 
another, as t crosses the value t*. Specifically, if the earlier solution 0„(t*, r) belongs to S^, 
then the later solution 6+(t*,r) belongs to S^^sn^- 

Equation ( |4.25|) shows that the change of the gauge field winding number SN^\t=t* causes 
the change of the spinor twist number at t = t* to be given by 



6K 



B:A\t=t* 



-6K 



(4.26) 



regardless of the existence of a fermion zero mode at t = t* (this contribution is labeled A). 

Before we continue with the evaluation of the change of spinor twist number, we need 
to address the continuity issue for the solution /i of Eq. ( [4.4| ). According to Eqs. ( [4.1 9|) 
and (|4.23 ), the time derivatives of the gauge-invariant X{t,r) and TZ{t,r) fields are not well- 
defined at the spacetime point {t*,r*). This implies that the function ji{t,r), as defined in 
Eq. ([4.5| ), does not have a well-defined value either. Note that all these problems can be 
traced to the discontinuity of the unitary matrix A in Eq. ( p.8|) , which, in turn, is caused by 
the ill-defined argument (p(t, r) of the field x(t, r) at the spacetime point (t*, r*) where x(t, ^) 
vanishes. By performing the inverse unitary transformation of Eq. ( p.7|) , \l/ — >^ = A\l/, 
one finds that the ji "expectation value" ( ^ 



is again given by 



OH 
'dt 



(t* 



(4.27) 



R 



with the implicit integral over r on the right-hand side running over [0, i?]. The right-hand 
side of Eq. ( [4.27] ) has a well-defined value at t = t* for smooth background fields ao(t,r), 
ai{t,r), and x{'t,r). This shows that the discontinuity of the function ji{t,r) at {t*,r*), 
caused by an ill-defined function (p{t, r), can be absorbed into the local change of the spinor 
twist number 6NQ-A\t=t* via the relation ( [4.26| ). 

With the discontinuity oi ji{t,r) at {t*,r*) absorbed into the local change of the spinor 
twist number 6NQ.j^\t=t*, the solution /i(t,r) of Eq. ([1.4|), explicitly given by Eq. ([4. 71), 
takes a well-defined value at t = t* and may produce a local change of spinor twist number 
that is not associated with the local change of gauge field winding number. 

Finally, we are ready to consider the additional effect on the local change of the spinor 
twist number due to the presence of a fermion zero mode at t = t* (this contribution will 
be labeled B). Since Jiit* ^ oo) as given by Eq. ( [4.27| ) has a well-defined value at t = t*, we 
have a unique time gradient for the level crossing. 



t=t*+e 



dE 
'dt 



Ji(r,oo). 



(4.28) 



t=t* 



Now, generic background gauge fields with j7i(t*,oo) ^ produce the following change of 
spinor twist number at t = t* [see Eqs. ( [4.14|) above]: 
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0]B\t=t* 



sgn [Ji(r,oo)] 



(4.29) 



in addition to the contribution SNQ-A\t=t* given by Eq. ( |4.26|) . [Remark that SNq-b = 
± 1 for the generic case considered.] According to Eq. ( [4.28|) , the level crossing at t = t* 
is determined by sgn [J7i(t*, oo)]. We therefore deduce the following relation between level 
crossing and the change of the spinor twist number: 



sgn 



'dE 




_ dF 


t=t*- 



6NQ\t=t* - 6N0.A\t=f 



(4.30) 



where the total change of the spinor twist number at t = t* is given by the sum of both 
contributions 



6K 



e\t=t* 



e;A\t=t* 



6K 



e:B\t=t* 



(4.31) 



Combining Eqs. ( [4. 261 ) and ( [4. 301 ), we find that the local spectral flow + e, t* — e" 
is given in terms of the local winding factor and twist factor, 



^[r + e, r - e] = 5N^\t=t* + 5Ne\t=t' 



(4.32) 



which is the main result of the present subsection. Here, and 6Nq are defined by Eqs. 
( [4.22| ) and (|4.13|) , respectively, and e is a positive infinitesimal. 



4- Special and generic gauge field backgrounds 



Let us, finally, discuss the case of having a fermion zero mode at t = t*, for which the 
"expectation value" J'i(t*,oo) vanishes, 

POO 

Ji(r,oo)= / dr |^(r,r)^9^7^-9^A sin20]|,^,, 
Jo 

= 0. (4.33) 

This implies dE/dt\t=t* = 0, according to Eq. ( [4.28| ). A nonvanishing d'^E/dt'^\t=t*, with 
vanishing first-order derivative, now corresponds to the absence of level crossing (the fermion 
energy eigenvalue just touches E = at t = t*). But a nonvanishing d^E/dt^\t=t*, with 
vanishing first- and second-order derivatives, again has level crossing. One therefore needs 
to find the first nonvanishing derivative of the energy corresponding to the fermion zero 
mode, in order to determine whether or not level crossing occurs. 

So far, we have considered only the solution \l/(t, r) of the zero-eigenvalue equation ( p.2aD , 
not the full spectrum ( p.la| ) of the time- dependent Dirac Hamiltonian. This makes it im- 
possible to obtain the relation between the local change of the spinor twist number and the 
time derivatives of E{t) att = t* beyond the leading-order approximation. Still, the relation 
( |4.32| ) can be shown to hold for generic SU{2) gauge field backgrounds, since background 
fields with J7i(t*, oo) = form a class of measure zero. 

Start from the gauge-invariant background fields X{t, r) and TZ{t, r) at the time slice t = t* 
where the fermion zero mode resides, with profile functions \E'(t*,r) and Q{t*,r). Under an 
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infinitesimal time shift t = t* ^ t* + e, generic X{t,r) and TZ{t,r) vary according to Eq. 
( |4.1|) , witli independent first-order coefficients (possible differences above and below t = t* 
are not important for the present argument). The change with time of these background 
fields is to first order in e 

6X = edtX\t=t*, 6n = edtTZ\t=t*. (4.34) 
Next, define the following "expectation value" of 6\ : 



<S\>t=t* = / dr|^(r,r)|^ SX{t,r)\,^,, 
Jo 

xsin2e(r,r), (4.35) 

with an integration measure weighted by the known function sin20(t*,r). The analogous 
"expectation value" of 671 is 

POO 

<6n>t=t*= dr |^(r,r)|2 57^(^,r)|,^,., (4.36) 
Jo 

but without extra weight function. 

Now, recall that the functional J'i{t*,oo) as given by Eq. (|4.33| ) is proportional to the 
difference of these "expectation values" in leading order, 

e Ji(r, oo)= <5n >t=t' - <5\ >t=t* 

+ 0(6^) . (4.37) 

Near the origin of the two-dimensional space spanned by X = < 5TZ >t=t* and Y = < 
SX >t=t*, the class of background gauge fields with j7i(t*,oo) = therefore coincides with 
the one-dimensional subspace 

{{X,Y)\X = Y}, (4.38) 

which is of measure zero. This shows that the quantity j7i(t*,oo) is nonzero for generic 
background gauge fields away from t = t* and that the relation between the local change 
of the spinor twist number and level crossing as given by Eqs. (|4.17 ) and (|4.3CI|) holds in 
general. 



B. Relation between spectral flow and SU{2) gauge field background 

The results of the previous subsection can be summarized as follows. For generic regular 
bosonic fields of the effective (l+l)-dimensional theory (|2.8| ), the spectral flow for 
the time interval [ti,tf] is given by the sum of the winding factor ( p.l7| ) and twist factor 

J^[tf,U] = AN^[tf,U] + ANe[tf,U]. (4.39) 
This result is simply the grand total of all level crossings ( [4.321) . 
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As mentioned before, the spectral flow T in Eq. ( [4.39| ) is defined as the number of 
eigenvalues of the effective Dirac Hamiltonian ( 3.1b|) that cross zero from below minus the 
number of eigenvalues that cross zero from above, for the time interval with ti < tf. 

The quantity is an integer, by definition. But the winding factor AA'^^ and the twist 
factor ANq also take integer values in general. 

Let us, nevertheless, discuss the special cases for which relation ([4.39|) is not applicable. 
The spectral fiow J-' from Eq. ( [4. 391 ) would not have a well-defined integer value if the gauge 
field winding number N^{t) or spinor twist number iVe(t) were ill-defined or noninteger at 
time slice t e Tij = {ti,tf}. In order to simplify the discussion, we exclude static field 
configurations from our considerations. 

The gauge field winding number N^{t), in particular, is not well-defined for a time slice 
t = t^^^ G Tif if the function x(t*^^\r) has a zero, see Eqs. ( |2.15| ) and ( |2.16| ). Assume 
that the zero of x occurs for t'^^^ = ti and that this is the only problem. In this case, one 
can simply choose a real number 6tf'\ so that the field xi^i''^) has no zero at the time slice 
t = tf'^ = U + Stf'K For the new time interval [tf\tf], one then obtains an integer- valued 
winding factor AN^[tf,t[^^]. The other case of having the zero of X at t(^) = tf can be 
treated in the same way. 

Alternatively, the spinor twist number N^lt) can take a half-odd- integer value for time 
slice t = t^"^^ G Tif. Now recall that a half-odd-integer spinor twist number implies the 
existence of a fermion zero mode; see the paragraph below Eq. (|3.19b ). This makes it 
impossible to properly define the spectral fiow for the exact time interval [ti,tf]. Assume 
that the zero mode occurs for t^"^^ = ti and that this is the only problem. In this case, one 
can choose a real number 6tl , so that A^e(^) takes a well-defined integer value at the time 
slice t = tf^ = ti + 6tf'^ . This is always possible, because a fermion zero mode corresponds 
to an asymptotically unstable solution 6 of Eq. ( |3.13a| ). For the new time interval [tf'\tf], 
one then obtains an integer- valued twist factor ANQ[tf,tf'^ ]. The other case of having the 
zero mode at t*-^-* = tf can be treated in the same way. 

Henceforth, we assume generic spherically symmetric SU{2) gauge field backgrounds, 
so that the right-hand side of Eq. (|4.39| ) is well-defined and the sum of two integers. It 
is, furthermore, clear that the winding factor AA^j^[tj,tj] is entirely determined by the 
background gauge fields [see Eqs. ( |2.15|) and ( |2.16|) ]. But also the twist factor AN^ltf.ti] 
can be expressed solely in terms of background gauge fields [see Eq. ( |3.22| )]. The relation 
( [4.39| ) thus connects a property of the fermions, the spectral flow JF, to two characteristics 
of a generic spherically symmetric SU{2) gauge fleld background, the winding and twist 
factors. This is the main result of the present paper. 



V. SPECTRAL FLOW FOR LUSCHER-SCHECHTER SU{2) GAUGE FIELDS 

In this section, we discuss the existence of fermion zero modes and the corresponding 
spectral flow for certain explicitly known time-dependent spherically symmetric solutions of 
the SU{2) Yang-Mills equations. This allows for a nontrivial check of the relation ( [4.39[ ) 
found in Sec. |VD 3| . Throughout this section and in the flgures, the same (arbitrary) mass 
scale is used to make the spacetime coordinates and energy dimensionless. 
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A. Brief review of the LS solutions 



The solutions considered in this section are spherically symmetric solutions of the SU (2) 
gauge field equations, which describe collapsing and re-expanding shells of energy. The 
corresponding (l+l)-dimensional field equations from the reduced action (|2.8a| ) read 



[-D' + i\xf-l)/r']x = 0. 



(5.1a) 
(5.1b) 



Remarkably, Liischer and Schechter were able to obtain analytic solutions of these coupled 
partial differential equations [ll2| , p!3| . 

The Liischer-Schechter (LS) solutions can be represented as follows (see Refs. [0,|1O| and 
references therein): 



a = Re X = 1 + Q'(t) cos^ w, 
P = lmx = (1/2) ^(t) sin2w. 



with the new coordinates 

r = sgn(t) arccos 



w = arctan 



1 - + 
2r 



(5.2a) 
(5.2b) 
(5.2c) 



(5.3a) 



(5.3b) 



Using the Ansatz (^.2|) , the field equations (|5.1| ) are reduced to a single nonlinear second- 
order differential equation for ^(t). 



d^g 
dr2 



+ 2g(g + l)(g + 2) = 0. 



(5.4) 



The ordinary differential equation (5^) can be interpreted as belonging to a mechanical 
system consisting of a particle trapped in a double-well potential V{q) = |g^(g + 2)^. The 
conserved total energy e of the particle trapped in the potential V is then 

2 



dq 
d7 



+ Viq), V{q)^q\q + 2f/2. 



(5.5) 



The general solution of Eq. ( p.4| ) depends on the energy parameter e and the r-translation 
parameter tq , together with a further discrete parameter = ± 1 . The solutions of Eq. ( ^.4|) 
can be divided into two classes, one with energy e < 1/2 and the other with energy e > 1/2. 



Explicitly, the LS solutions are ||T^ , p!3 

g(r) = -1 + C (1 + y2i)V2dn [(1 + y2i)i/2(r - Tq) | 



for e < 1/2, 

g(r) = -1 +((1 + V2i)i/2cn [(8e)i/4(^_^^)|^j ^ > 1/3, 



(5.6a) 
(5.6b) 
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with the modulus defined by 



m 



2V2e 



(5.7) 



Here, dn['u|m] and cn[M|m] are Jacobi elhptic functions [^PH . 

For e < 1/2, there exists no spacetime point where xit^f) vanishes. For e > 1/2, on the 



other hand, there are zeros of x{i}^) [|T0 | 

l + 2n 



tn = tan To + 



Km] 



(8e)V4 

where n is an integer that satisfies the condition 



TT 

- < 

2 - 



To 



l + 2n 
(8^ 



K{m] 



vr 

< -. 

- 2 



(5.8) 



(5.9) 



Here, K[m) is the complete elliptic integral of the first kind. 



K{m) 



fdu [{I- 
Jo 



mu 



-1/2 



(5.10) 



The existence of a spacetime point {t*, r*) where x{t, r) vanishes is, in general, associated 
with the change of the winding number N^{t). It has been shown in Ref. [0 that the change 
of the winding number N^{t) plays an important role in fermion number violation (see also 
Ref. [|l^ for related results). Indeed, we have studied several LS solutions with e < 1/2 and 
found that the necessary condition ( p.25| ) for the existence of a fermion zero mode is never 
satisfied. In the following, we shall therefore only consider LS solutions with e > 1/2. But 
before we turn to the fermion zero modes, we mention one particular aspect of the LS gauge 
field background for e > 1/2. 



B. LS quasi-sphaleron 

For energy parameter e > 1/2, the field x(t, r) has at least one zero at a particular 
spacetime point. In order to simplify the analysis, we take for the ( and tq parameters in 
the solution (|5.6b| ) the following values: 



C = +l, ro = -{8e)-'/'K{m), (5.11) 

with m defined by Eq. ( |5.7| ). For the choice of Tq from Eq. ( p.llj ), one of the zeros of x{t, f) 
occurs at the time slice t = to = 0. Also note that for this tq the LS solution has time-reversal 
(anti)symmetry, namely p{—t,r) = p{t,r) and 7l{—t,r) = —7l{t,r), with p{t,r) and 7l(t,r) 
defined in Eq. (|3.14| ). In the following, we consider the t = time slice of these particular 
LS background gauge fields [parameters ( = +1, To from Eq. ( |5.11| ), and e > 1/2], for which 
the zero of x(0,t) occurs at r = 1, according to Eq. (|5l8| ). 
The LS gauge field background at t = is represented by 

a^(0,r) = 0, x(0,r) = sinu;(0,r), (5.12) 
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up to a U{1) gauge transformation Given that the real function x(0,t) vanishes and 

changes sign at r = 1, the configuration of Eq. ( ^.12 ) quahtatively resembles the sphaleron 
solution of the electroweak SU{2) Yang-Mills-Higgs theory For this reason, we call 

the configuration given by Eq. (|5.12| ) the "LS quasi-sphaleron." Note that the LS quasi- 
sphaleron does not satisfy the static field equations, since the energy changes under a scale 
transformation of the fields (there is no natural mass scale for classical Yang-Mills theory). 

We will now show that this LS quasi-sphaleron corresponds to the top of a potential en- 
ergy barrier which separates configurations with different A^^, just like the electroweak spha- 
leron 1^. The energy functional [[7|,pT![ for (H-l)-dimensional gauge field solutions can, in 
fact, be written as 



E = EKit) + Ep{t), 



with 



EK{t) 



9 Jo 



TC 



Ep{t) = -r 



9 Jo 



dr 



dr 



(5.13) 

(5.14a) 
(5.14b) 



where p{t,r) equals | x(^)'")l and ip{t,r) is the (bosonic) gauge-invariant field defined by 

-2e^,^(t,r) = rV^,(t,r). (5.15) 

Here, we have divided the energy into the kinetic part E^ and potential part Ep. The 
reason for putting the ■^/'-dependent terms into the kinetic part of the energy is that, for the 
gauge choice oq = 0, the scalar field ip becomes proportional to the time derivative of the ai 
field, namely = — r^^octi. 

Consider then the potential energy Ep{t) associated with our particular LS gauge field 
solutions. Using the basic properties of the Jacobi elliptic function cn[M|m], one can 
prove that the potential energy Ep{t) has a local maximum at t = 0, that is, for the LS 
quasi-sphaleron configuration. Furthermore, the conserved total energy of the LS solution 
with To parameter ( ^.11| ) is given by 



E = 2eE. 



quasi— sph) 



(5.16) 



in terms of the energy parameter e and the static energy of the LS quasi-sphaleron, -Equasi-sph 
= Ep{0). 

Figure |1] shows the time development of Ep{t) for the LS background with e = 1, ( = +1, 
and To ~ —1.4271 from Eq. (|5.11| ). The corresponding topological charge Q ~ —0.70 is 
noninteger, see Ref. ^ for further details. More importantly, the potential barrier of Fig. |I| 
separates two regions with different winding number (AA^^ = —1). The LS quasi-sphaleron 
at t = resembles in this respect also the sphaleron of the electroweak Standard Model . 

For comparison, consider the LS gauge field background given by the trivial solution 
q = —1 of Eq. (|5.4|) , with the (l-(-l)-dimensional field configurations 



ip{t,r)=0, x(t, r) = sin ti?(t, r). 



(5.17) 
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as follows from Eqs. (|2.11| ), ( p.2|) , and ( p.l5|) . This corresponds to the de Alfaro, Fubini 
and Furlan (AFF) solution [^. Note that the AFF gauge field coincides with the LS quasi- 
sphaleron ( |5.12| ) at t = 0. On the other hand, the AFF solution has complete time-reversal 
symmetry and the kinetic energy Exit) as given by Eq. (|5.14aj) is zero at t = (see Fig. 
0). This result suggests that the AFF solution provides the time-dependent gauge field 
solution with minimum total energy to form the LS quasi- sphaleron. In other words, the 
AFF gauge field simulates an "imploding and exploding LS quasi-sphaleron;" cf. Ref. P3 



(For the electroweak sphaleron, the dynamics of the gauge and Higgs fields has been studied 
numerically. See, for example, Refs. 



C. Fermion zero mode of the LS quasi-sphaleron 



We now turn to the fermion zero-eigenvalue equation ( |3.9| ) for the t = LS gauge field 
of the previous subsection, i.e. the LS quasi-sphaleron (|5.12|) . Note that the gauge-invariant 
function A(0,r) = p{0,r)/r is nondifferentiable at r = 1, which is the only point where 
p(0,r) = |x(0,t)| vanishes. For this reason, we introduce a differentiable field A(0,r), 
defined by 

~X{0,r) = K{0,r)/r, (5.18a) 
K(0,r) = sgn(l - r) |x(0,r)|. (5.18b) 

Then, x(0,r) can be represented by 

X(0,r) = «:(0,r) exp[^ (^(0, r) ], (5.19) 

where (p{0, r) is a differentiable function of r. 

It is a simple exercise to verify that the LS field k(0, r) has the following inversion 
symmetry: 

K{0,l/r) = -K{0,r). (5.20) 

This inversion symmetry, most likely, traces back to the conformal symmetry transformation 
x'^ — i> x^/x"^ of classical Yang-Mills theory. Without loss of generality, we consider, in the 
following, smooth background fields with k{0, 0) = 1 and (^(0, 0) = 0. 

In terms of the differentiable fields A and TZ = (ai — 9r<^)/2, one finds that the zero-energy 
fermion equations ( p.l3| ) at t = become 

drQ = -A sin 29 + TZ, (5.21a) 

= A|^| cos2e, (5.21b) 

with boundary conditions 

9(0,0) = 1^(0,0)1 = 0. (5.22) 

Using the inversion symmetry ( [5. 20] ) and the asymptotics ( |3.11D of the rotator, it follows 
from Eq. ( ^.21b| ) that there exists a fermion zero mode if and only if limr^oo ©(0, r) = Nn, 
with an integer. 
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We can now explicitly construct the fermion zero mode for the LS background at t 
Using the field equations ||7,10|, it is relatively straightforward to show that 



0. 



n{0,r) oc 



dr2 



oc cn[i^(m) 



m 



0, 



(5.23) 



where the last identity can be found, for example, in Ref. [Note that vanishing TZ 



does not contradict the necessary condition (|3.25|) that was derived for Eqs. ( |3.13|) with a 
function A > 0, whereas Eqs. ( [5. 211 ) have a function A that changes sign.] 
For TZ = 0, the solutions to Eqs. ( p.21j ) are simply given by 



e(o,r) = o, 

1^(0, r)| = 1^(0, ro)| exp 



dr' A(0,r') 



(5.24a) 
(5.24b) 



The argument following Eq. ( [5.22| ) ensures that the solution represented by Eqs. ( p.24| ab) 
is normalizable, which completes the construction of the fermion zero mode. The inversion 
symmetry ( ^.20|) , together with the result (|5.23| ), provides a sufficient condition for the 
existence of the fermion zero mode. 

For e > 1/2, the fermion zero mode amplitude |\l/(0,r)| of Eq. (|5.24b|) is shown in Fig. 
^, with an arbitrary normalization. Specifically, the (3 + l)-dimensional fermion zero mode 
is purely left-handed and given by Eqs. (|2.6| ) and (|2.9| ), with the two-component spinor 

= \E' from Eqs. ( |3.7| ), ( p.l2| ), and ( p.24| ). These last equations are to be evaluated with 
the functions K(0,r) and <^(0,r) defined in Eq. (|5.19| ). [Of course, the solution can also be 
obtained from Eqs. ( ^.13| ), which are given in terms of A(0,r) and (p{0,r). In this case, the 
function 0(0, r) is found to have a step function at r = 1, but so does the transformation 
matrix A from Eq. (|3.8|) .] 

To summarize, the LS quasi- sphaleron at t = has a chiral fermion zero mode and 
resembles in this respect the electroweak sphaleron which also has a chiral fermion zero 
mode (see Refs. [p6| -p8[1 and references therein). The fermion zero modes of the LS quasi- 
sphaleron and the electroweak sphaleron are qualitatively the same. Moreover, there is 



spectral fiow associated with both the electroweak sphaleron (see Refs. [P7|-[2m) and the LS 
quasi-sphaleron (see Sec. V D 3 below). This behavior differs from that of the AFF solution 
(an "imploding and exploding LS quasi-sphaleron"), for which the fermion zero mode exists 



at all times [23|, without level crossing. 



D. Level crossings for large energy parameter e 

Throughout this subsection, we consider the specific LS gauge field solution with pa- 
rameters e = 20, ( = +1, and tq ^ —0.54197 from Eq. ( 5.11 ). Figures ^ and ^ give the 



behavior of the potential energy Ep{t) of this solution. The corresponding topological charge 
Q —0.13 is noninteger; cf. Ref. 0]. The LS quasi-sphaleron at t = is the same as for 
the e = 1 case, but no longer corresponds to a global maximum of Ep{t). 

In order to determine the spectral fiow for this particular gauge field background, all 
fermion zero modes need to be determined, to which we turn first. 
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1. Extra fermion zero modes from changes in winding number 



According to Eq. ( ^.8|) , there are zeros of this LS field x{'t,f) at the spacetime points 
(t-i,r-i) ^ (-1.889,2.137), (to,ro) = (0,1), 

(i+i, r+i) ^ (+1.889, 2.137). (5.25) 
The winding number N^{t) takes the following values in the different time regions: 




for t G (— oo, t_ 



1 



^ ^ -1 for t G (to,t+i), ^ ^ 

for t G (t+i, +oo). 

This gives the global winding factor 

AAr^[+oo,-oo] = - (-1) = 1. (5.27) 

As discussed in Sec. V C, there exists a fermion zero mode a.t t = to = (see Fig. 
But there are two more fermion zero modes precisely at t = t_i and t = t^i. Three 
preliminary steps are necessary for the proof. 

First, define smooth fields X(tn,r) and lZ{tn,r) at the time slices t = tn, for n = ±1, 

~X{tn,r) = K{tn,r)/r , (5.28a) 
Tlitn^r) = [ai(t„,r) -9,^(t„,r)]/2 , (5.28b) 

with the differentiable function n(tn,r) = p(tn,r) x sgn(r„ — r) and the smooth argument 
(p{tn,r) of the Higgs-like field x = exp[2(^]. 

Second, perform a scale transformation x = r/rn-, so that the radial point r = r„ where 
the field xi'^m'^) vanishes corresponds to a; = 1. 

Third, establish that the Liischer-Schechter solution for the parameters chosen has the 
following inversion symmetry at fixed time slices t = tn'- 

A(t„, 1/x) = -x^ X{tn, x) , (5.29a) 
n{tn, 1 /x) = -x^ n{tn, x) , (5.29b) 

for n = ±1. The existence of this inversion symmetry has been verified analytically with 
the help of Mathematica 4.0 p]] . 

After these preliminaries, we turn to the possible existence of fermion zero modes at the 
time slices t = t±i. Consider the zero-energy fermion equations given by Eqs. ( p.l3| ) at 
t = tn, for n = ±1, with the smooth background fields A(t„,r) and 7?.(t„,r). The chiral 
Yang-Mills theory ( |2.1| ) is scale invariant and the rescaling x = r/r±i at t = t±i does not 
alter the structure of Eqs. ( p.l3|) : 

^ = -A(x) sin 2e(x) + n{x) , 6(0) = , (5.30a) 
dx 



d|^| 
dx 



A(x) |^(x)| cos20(x) , 1^(0)1 = 0, (5.30b) 
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with the dependence on t±i temporarily dropped. The differential equations ( p. 301) are 
symmetric under the inversion transformation x ^ 1/x and so are their solutions 0(x) and 
The inversion symmetry implies 

lim e = lim e = , lim |^| = lim |^| = , (5.31) 

x^O X— >oo x^O x^oo 

with 1^1 oc 1/x for large x. This shows that there exist fermion zero modes at both t = t^i 
and t = t+i. 

The inversion symmetry ( |5.29| ), after the appropriate scale transformation, provides again 



a sufficient condition for the existence of fermion zero modes at t = t±i. For Liischer-Schech- 
ter gauge field backgrounds with arbitrary tq and e > 1/2, the inversion symmetry ( [5.29[ ) 



holds, in fact, at any time slice t = where xi't, r) has a zero. This then proves the existence 
of fermion zero modes at all t„. 

Figure ^ gives the profile functions of the fermion zero mode at t = t_i, obtained from 
the numerical solution of Eqs. ( ^.30| ). The profile functions of the fermion zero mode at 
t = t+i are identical, except for a change of sign of 0. These functions are quantitatively 
different from those of the LS quasi-sphaleron [see Eqs. (|5.24|) and Fig. ^, but qualitatively 
the same. 



2. Extra fermion zero modes from changes in twist number 



In order to locate all possible level crossings, we are guided by the change of the spinor 
twist number NQ{t). The spinor twist number NQ{t) takes the following values in the 
different time regions: 



+1 for t e (-00, -ta), 
for te{-t,,+ta), 
-1 for te(+ta,+oo), 



Ne{t) = { for tG (-t„+t„), (5.32) 
with the numerical estimate ~ 2.924. The corresponding twist factor is thus given by 



AiVe[+oo,-oo] = -1 - 1 = -2. 



(5.33) 



In addition to the fermion zero modes at to aiid t±i, which are associated with the 
change of the gauge field winding number N^{t), there exist two more fermion zero modes 
precisely at t = dita- The analysis of these fermion zero modes is straightforward and 
the normalizability condition is found to hold, provided Q{±ta,r) approaches a half-odd- 
integer multiple of TT as r — i> oo; see Eq. ( p.l3b|) . Note that the exact value of ta is defined 
implicitly by the relation lim/j^oo ©(i^a^-R) = T7r/2, where the solution Q{±ta,R) of the 
differential equation ( |3.13a] ) with boundary condition ( p.l5a| ) can be obtained by the method 
of successive approximations |[l^J15[| . 

Figure ^ gives the profile functions of the fermion zero mode at t = —ta, obtained from 
the numerical solution of Eqs. ( 3.13 ) and ( 3.15| ). The profile functions of the fermion zero 
mode at t = +ta are identical, except for a change of sign of 6. Figures |^ and ^ show the 
time-variation of the solutions Q{t,r) of the transformed Riccati equation ( p.l3a| ) around 
t = ±ta, which demonstrates that the fermion zero modes at t = ±ta sit at bifurcation points 
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for different N^^s. [These results provide an example for the general discussion leading up 
to Eqs. (|4.14| ) in Sec. IV A 2.] Obviously, Figs. || and || are related, because of the 
time-reflection properties of the background fields mentioned below Eq. ( |5.11| ). 

The fermion zero modes at t = ±ta are qualitatively different from the ones at to and t±i 
(compare Fig. |^ with Figs. ^ and ^). These fermion zero modes occur, in fact, for Higgs- 
like fields x(±ta,r) without zeros. This differs from the cases discussed in the literature 
T(]| , P^J5(]| . Apparently, the long-range behavior of the background SU{2) gauge fields plays 



a crucial role for the existence of these extra fermion zero modes (see also the discussion in 
Sec. VI). 



3. Spectral flow 



For the fermion zero modes at t = — t^, t_i, t^, t+i, and +ta, we calculate the following 
time-gradients of the energy eigenvalue of the effective Dirac Hamiltonian: 



dE 


dE 


'dt 


t=-u t 


dE 


dE 


dF 


t=t., , 



t=+ta 



t=t4 



-0.03 < 0, 



+0.08 > 0, 



dE 
dF 



-5.00 < 0. 



(5.34) 



t=to 



In addition, we have checked that there are no further fermion zero modes, at least for the 
time interval [-200, +200]. 

The level crossings corresponding to Eq. ( |5.34| ) give the following value for the spectral 
flow (starting from t = —ta and ending at t = +ta): 



jF[+oo, — oo 



1 + 1 - 1 = -1. 



(5.35) 



It would certainly be interesting to calculate the spectrum of the time-dependent effective 
Dirac Hamiltonian numerically (cf. Ref. [^]) in order to see whether or not the pattern 
( ^.35| ) corresponds to a single energy level crossing E = five times. Anyway, the total 
value ( |5.35 ) agrees with the spectral flow obtained from the relation ( |4.39|) : 



jF[+oo, — oo 



AiV- 
1 - 



-oo, — oo ] + AA^0[+oo, — oo 
= -1, 



(5.36) 



where the results ( ^.27| ) and ( |5.33| ) have been used. This demonstrates the role of the twist 
factor ANq for the spectral flow, which has not been noticed before to our knowledge. 

Recall that the explicit results of the present subsection are for the particular LS gauge 
field background with energy parameter e = 20, together with ( and tq from Eq. ( p. 11 ). 
The spectral flow for LS background gauge fields turns out to be solely given by the winding 
factor AA^^ if e is smaller than e* ~ 5.37071, which is the numerical solution of the following 
equation: 
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(5.37) 



Briefly, the argument runs as follows. First, the numerical results (and a heuristic argument) 
give ta > i+i- Second, the times t±i move toward ±00 as e approaches e* from above; see 
Eqs. ( |5.8| ) and ( ^.11| ). Together, this implies that the contribution of the twist factor to 
the spectral flow vanishes for e < e*, giving jF[+oo, —00] = AA^^[+oo, —00]. (The winding 
factor AA^^[+oo, —00] is, of course, identically zero for e < 1/2.) 

To summarize, the twist factor AA^^e can play a significant role for the spectral flow in 
certain LS gauge field backgrounds, provided the energy parameter is large enough (e > e*). 



VI. DISCUSSION 



In this paper, we have studied real-time anomalous fermion number violation by di- 
rectly investigating the zero-eigenvalue equation ( |3.1a|) of the time- dependent effective Dirac 



Hamiltonian for spherically symmetric massless chiral fermions and SU (2) Yang-Mills gauge 
field backgrounds. For these spherically symmetric SU{2) gauge field backgrounds, we 
have found a relation between the spectral flow and two characteristics of the gauge fields. 
Physics applications of this result are based on the assumption that anomalous production 
of fermions is confined to the spherically symmetric partial wave; cf. Ref. Perhaps the 



most important application would be for electroweak baryon number violation in the early 
universe, in particular at temperatures above the electroweak phase transition (see Ref. fSlI 
for a review). 



Since we adopt an approach different from the one of previous work |5|,|9|Jl0[| , we are able 



to observe certain new features of real-time fermion number violation within the spherically 
symmetric Ansatz. These features include the spinor twist number NQ{t) obtained from the 
Riccati equation ( |3.16| ) at a single fixed time t and the corresponding twist factor AN0[tf, ti ] 
which is the change of the spinor twist number over the time interval [ti,tf], with ti < tf. 
Relation ( [4.39| ) then gives the spectral flow JF as the sum of this twist factor AA'^e and the 
winding factor AA^^, which is the change of the Chern-Simons number of the associated 
vacuum sectors of the gauge field background. Mathematically, the relation ( [4.39| ) takes the 
form of an index theorem, restricted to spherically symmetric fields; cf. Ref. 0]. 

In order to get better insight into the meaning of this relation ( |4.39| ), we have investigated 
level crossings for a particular class of Liischer-Schechter (LS) gauge field solutions |l2| , p!3 



The results clearly demonstrate the role of the twist factor for the spectral flow. See, in 
particular, Eqs. ( [5.35| ) and (|5.36| ) . The nonvanishing global effect of the twist factor on the 
spectral flow is partly due to the fact that the LS gauge fields form propagating solitons in 
the effective (H-l)-dimensional theory; cf. Refs. [[7|,pi|. The fields are thus nondissipative in 



the (l-l-l)-dimensional world. For such background gauge fields, the field-theoretic approach 
adopted in Refs. P,pi!| is, strictly speaking, not applicable. 



Although nondissipative in the (l+l)-dimensional context, these LS gauge field solu- 
tions are dissipative in (3+l)-dimensional spacetime. The (3+l)-dimensional energy density, 
which is obtained from the (l+l)-dimensional energy density ( |5.14| ) by dividing by 47rr^, 
approaches zero uniformly for early and late times (t ±00). 
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At this moment, we propose to classify dissipative spherically symmetric SU{2) gauge 
field solutions into two categories. A spherically symmetric SU{2) gauge field solution 
is called strongly dissipative, if both the (3+l)-dimensional and (l+l)-dimensional energy 
densities approach zero uniformly for large times (t — > ±00). On the other hand, a spherically 
symmetric SU{2) gauge field solution is called weakly dissipative, if the (3+l)-dimensional 
energy density dissipates with time, but not the (l+l)-dimensional energy density. Note 
that the LS gauge field solutions considered in Sec. V D are weakly dissipative, according 
to this terminology. 

For strongly dissipative spherically symmetric SU{2) gauge field solutions, the rotator 
TZ{t, r) in the transformed Riccati equation ( |3.13a| ) lacks the strength to give a nonzero 
spinor twist number A^^e at large times (see Appendix A for the proof). The relation (f4.39|) 



then predicts that the spectral fiow is solely given by the winding factor AA^^, which 
reproduces the known result 0,0. An isolated change of spinor twist number can still 
contribute to the local pattern of level crossing ( |4.32| ) . 

For weakly dissipative or nondissipative spherically symmetric SU{2) gauge field solu- 
tions, a nonvanishing twist factor AA'e can even make a global contribution to the spectral 
flow JF. As mentioned above, this has been verified for certain LS solutions [jl2|,|13|. This 
behavior does not follow directly from the perturbative triangle anomaly , which detects 
only the (noninteger) topological charge. Recall that the standard perturbative calculations 
(Feynman diagrams) essentially neglect the interactions of incoming and outgoing particles. 



i.e. the interactions are "turned off" in the asymptotic regions 



To summarize, the spectral flow result ([4.39|) does not assume strongly dissipative spher- 



ically symmetric SU{2) Yang-Mills gauge fleld backgrounds, in contrast to previous studies 
||5|j9|JT0[| . It may, therefore, be applied to real-time anomalous fermion number violation for 
weakly dissipative or nondissipative spherically symmetric SU{2) Yang-Mills gauge fleld 
backgrounds. Moreover, preliminary results indicate that the relation (|4.39|) can be adapted 
to the case of chiral fermions interacting with spherically symmetric SU (2) Yang-Mills and 
Higgs flelds. 

The main outstanding problem is, of course, to understand the role of the (appropri- 
ately generalized) twist factor in the full (3+l)-dimensional SU{2) Yang-Mills theory, not 
just the subspace of spherically symmetric conflgurations. Also, the corresponding index 
theorem needs to be established for the long-range Yang-Mills gauge flelds considered; cf. 
Ref. and references therein. Finally, the proper deflnition (if at all possible) of the 



second-quantized fermion number operator in general nondissipative Yang-Mills gauge fleld 
backgrounds requires further study. 



APPENDIX A: SPINOR TWIST NUMBER AND STRONGLY DISSIPATIVE 

SU{2) GAUGE FIELDS 

In this appendix, we calculate the asymptotic spinor twist number A'e(t), deflned by Eqs. 
( p.l3a| ), ( p.l5a| ), and ( p.20| ) in the main text, for strongly dissipative spherically symmetric 
SU (2) gauge fleld solutions. 

From the (l-l-l)-dimensional energy density ( 5.14|) and the corresponding fleld equations 



(see, in particular Eqs. (3.8) and (3.10) of Ref. [[T^]), the condition of strong dissipation 
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implies that there exists a small positive quantity e(t) <^ 1 at large \t\ , so that 

|r'^ - A(t,r)| < e(t) and |7^(^, r)| < e(t) , (Al) 
for < r < oo, together with the limit 

lim e(t) = . (A2) 

\t\~*oo 

See Eqs. ( p.9|) , ( p.l5|) , and ( p. 141) in the main text for the definition of X{t,r) and 7l{t,r). 
From the bounds ([Al|), we immediately obtain 



< r"^ - 2 e(t) < A(t, r) - |7^(^, r) | 



(A3) 



for < r < _Re = [2e(t)]^^. Since regular finite-energy gauge fields obey lim,.^oo'7^/A = 
[see Eq. ( p. 11 ) in the main text] , we have from Eq. ( |A3|) the following inequality for arbitrary 
r: 



0< \n{t,r)\<\{t,r), (A4) 

provided |t| is sufficiently large. 

With the inequality ( |A4D in hand, we are able to establish that the spinor twist number 
Neit) vanishes asymptotically. This can be shown by contradiction. Consider the trans- 
formed Riccati equation, 

dre{t,r) =V{t,r) +n{t,r), e(t,0) = 0, (A5a) 
V{t, r) = -A(t, r) sin 2e{t, r) , (A5b) 

at sufficiently large |t|. Now, assume that there exists a time slice t = t where the background 
fields obey inequality (|A4|) and that the solution 9(t, r) of Eq. (|A5|) belongs to the class 
SN{i), with positive integer (the case of negative A^ will be dealt with later). See Eq. 
( P.19a| ) in the main text for the definition of the solution class Sn. 

Strong dissipation gives |x(^, '^)| 7^ 0, so that the background fields X(t,r) and 7l(t,r) 
are smooth. This implies that the solution Q{i,r) G S]\f{t), with A^ > 1, is continuous and 
must cross the value +7r/4 at least once. Define r+ to be the largest radial distance for 
which Q{i,r) = + 71/4. In order for B(t, r) to reach the asymptotic value Ntt, the slope of 
r) at r = r+ clearly must be nonnegative, 

9^e(t,r+)>0, (A6) 

as long as 6 G SN(t), with A^ > 1. 

On the other hand, Eqs. ( [A4D and 
give the following inequality: 

drOit, r+) = - A(t, r+) + 7^(t, r+) < . (AT) 

This last result contradicts the earlier result ( [A6| ) , which was based on the assumption that 
0(t, r) G SN(t), with A^ > 1. Hence, Q(t,r) ^ S'Ar(t), for positive integer A^. The case of 



A5), together with the fact that sin29(t,r+) = +1, 
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negative integer N is ruled out in the same way. The conclusion is thus that Q{t, r) belongs 
to Soit). 

For strongly dissipative spherically symmetric SU{2) gauge field solutions, we find that 
the rotator TZ{t, r) at large times lacks the strength to overcome the resistance of the deviator 
T>[t, r), so as to give a nonzero spinor twist number at large times. In short, we have 

lim Ne{t) = lim A^e(t) = 0. (A8) 

This result shows that the twist factor AAi'el+oo, — oo] = Nq{+oo) — Nq{—oo) = does not 
contribute to the spectral fiow ( [4. 391 ), at least for the case of strongly dissipative spherically 
symmetric SU{2) gauge field solutions. 
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FIGURES 




FIG. 1. Time-development of the total energy E and potential energy Ep(t) for the Lii- 
scher-Schechter (LS) gauge field solution (0), (^.21 ), ( ^.6b ) with parameters e = 1, C = +lj 
and tq ~ —1.4271 from Eq. (Islll) . The potential energy Ep(t), given by Eq. (|5.14bD , is differen- 
tiable for all times and has a global maximum at t = 0. The configuration at t = is called the 
LS quasi-sphaleron. 
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FIG. 2. Time-development of the total energy E and potential energy Ep{t) for the de Alfaro, 
Fubini and Furlan (AFF) gauge field solution (U), (|^), (|5.17| ). The AFF configuration at t = 
coincides with the LS quasi-sphaleron at t = in Fig. 
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FIG. 3. Profile function |'I'(0, r)| of the fermion zero mode ( 5.24b| ) of the LS quasi-sphaleron, 
which corresponds to the t = configuration of Fig. ||. The dashed curve gives |^'(0,r)| with an 
arbitrary normahzation. The inversion symmetry r ^ 1/r is made manifest by use of the compact 
radial coordinate ry = (r — l)/(r + 1). 
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FIG. 4. Time-development of the total energy E and potential energy Ep[t) for the Lii- 
scher-Schechter gauge field solution (^^), ( |5.2| ), ( 5.6b| ) with parameters e = 20, C = +1, and 
To w -0.54197 from Eq. ( |5ll| ). See Fig. | for a close-up of Ep{t) near t = 0. 
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FIG. 5. Same as Fig. ^. Close-up of the potential energy Ep{t) near t = 0. The LS 
quasi-sphaleron at t = is only a local maximum of Ep{t), unlike the case of Fig. ||. 
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FIG. 6. Numerical solutions for the profile functions 0(t_i,r) and |^'(t_i,r)| of the fermion 



zero mode at t = t_i —1.889 for the Liischer-Schechter background gauge field (|2.4| ), (pi 
( |5.6b| ) with parameters e = 20, C = +1, and tq ~ —0.54197. The solid curve corresponds to 
G(t_i,r) X 2/7r and the dashed curve to |^'(t_i,r)| with an arbitrary normalization. The inversion 
symmetry x — > 1/x, with x = r/r_i, is made manifest by use of the compact radial coordinate 
r] = (x — l)/{x + 1). 
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FIG. 7. Same as Fig. ^, but for the fermion zero mode at t = —ta ~ —2.924. The soHd cm've 
corresponds to @{—ta, r) x2/tt and the dashed curve to |^'(— ta, r)| with an arbitrary normalization. 
Both functions are plotted against the compact radial coordinate r] = (j; — l)/(x + l), with x = r/2. 
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FIG. 8. Numerical solutions @{t, r) of the transformed Riccati equation ( 3.13a ), with boundary 
condition 0) = 0, at different times around t = —ta ~ —2.924 for the Liischer-Schechter back- 
ground gauge field (|2.4D, (|5.2D, ( p.Gb] ) with parameters e = 20, C = +li and tq ~ —0.54197. For 
t = —ta, there is a normalizable fermion zero mode, with limr__»oo ©(— ^a, = 7r/2 (see Fig. |^. 
The different solid curves for r > 10, from top to bottom, correspond to t + ta = —0.50, —0.05, 
—0.01, 0.00, 0.01, 0.05, and 0.50. With increasing time, the Q values at large r move toward the 
constant G value at r = 0. 
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FIG. 9. Same as Fig. ^, but for times around t = 2.924. For t = ta, there is a normalizable 
fermion zero mode, with hm^-^oo Qita,r) = — 7r/2. The different sohd curves for r > 10, from top 
to bottom, correspond to t-ta = -0.50, -0.05, -0.01, 0.00, 0.01, 0.05, and 0.50. With increasing 
time, the Q values at large r move away from the constant B value at r = 0. 
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